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SECTION  1 
INTRODUCTION 

1.1.  Purpose  of  the  Calculation  Series. 

This  report  describes  a  series  of  calculations  undertaken  to  compare  the  effects  of  geo¬ 
metric  approximations /simplifications  commonly  made  when  setting  up  hydrody¬ 
namic  code  simulations  of  blast  and  flow  in  tunnel  systems.  The  approximations  are 
made  to  save  computer  resources  by  running  the  calculation  in  two,  rather  than  three, 
dimensions.  A  two-dimensional  calculation  requires  less  computational  time  by  a  factor 
of  one  to  two  orders  of  magnitude.  It  also  requires  less  memory  to  run  and  less  storage 
space  for  the  computational  output.  In  addition,  setup  and  analysis  times  are  reduced. 
The  questions  to  be  answered  by  this  set  of  calculations,  which  included  two  different 
two-dimensional  representations  as  well  as  a  three-dimensional  representation,  are 
whether  the  approximations  provide  results  that  are  sufficiently  close  to  those  provided 
by  calculating  in  three  dimensions,  and  if  so,  which  two-dimensional  approximation  is 
the  best. 

The  three  calculational  representations  are  set  up  in  the  coordinate  systems  available  in 
SHARC  (S-Cubed  Hydrodynamic  Advanced  Research  Code),  a  state-of-the-art,  second- 
order-accurate,  multi-material  Eulerian  hydrocode.  These  include  a  two-dimensional 
cylindrical  system,  a  two-dimensional  Cartesian  (rectangular)  system,  and  a  three-di¬ 
mensional  Cartesian  system.  Each  representation  includes  an  explosive  in  a  detonation 
chamber,  an  input  tunnel,  an  expansion  chamber,  and  an  output  tunnel  connecting  to 
the  exterior.  For  each  chamber,  cross-sectional  areas  and  chamber  volumes  are  the  same 
in  each  configuration.  Yield  of  the  explosive  source  must  also  be  conserved  among  the 
different  representations.  How  this  was  accomplished  is  explained  in  detail  in  the  next 
section;  however,  the  two-dimensional  representations  are  not  unique,  and  some  alter¬ 
natives  are  discussed  in  Section  6.  All  of  the  calculations  assume  perfectly-reflecting, 
perfectly  smooth  tunnel  wall  surfaces.  Hence,  there  are  no  losses  to  the  walls. 
Gravitational  and  turbulence  effects  are  likewise  not  included. 

V 

1.2.  Geometric  Approximations. 

Three  calculational  configurations  were  specified,  described  as  two-dimensional 
Cartesian,  two-dimensional  axisymmetric  (cylindrical),  and  three-dimensional 
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Cartesian.  The  sketch  provided  by  the  Contract  Technical  Monitor  to  describe  the  con¬ 
figurations  is  included  as  Figure  1. 

The  detonation  chamber  was  given  as  a  cylindrical  chamber,  50.20  m  in  length  and  5.64 
m  in  diameter.  This  was  to  be  connected  to  the  input  tube  shown  in  Figure  1  for  running 
the  various  versions  of  the  problem.  Within  the  detonation  chamber,  a  cylinder  of  TNT, 
25.38  m  long  and  with  a  radius  of  1.41  m,  was  placed  coaxially  with  its  end  at  the  re¬ 
mote  end  (away  from  the  input  tube)  of  the  detonation  chamber.  The  explosive  was 
detonated  at  this  remote  end.  Using  a  density  of  1.56  gm/ cm3  for  TNT  with  the  calcu¬ 
lated  volume  of  the  explosive  cylinder,  158.52  m3,  the  amount  of  TNT  in  the  chamber 
was  247.29  metric  tonnes. 


Output 


Two-Dimensional  Cartesian  Two-Dimensional  Axlsymmetric 


Three-Dimensional  Cartesian 


Figure  1.  General  configurations  for  expansion  chamber  and  input  and  output  tunnels. 

Additional  information  provided  that  the  input  and  output  tunnels  were  30  m  in  length, 
8  m  wide,  and  6  m  high.  The  expansion  chamber  was  15  m  long,  70  m  wide  and  8  m 
high.  The  lateral  offset  between  the  input  and  output  in  the  Cartesian  configurations 

was  to  be  set  at  58  m. 
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1.2.1.  Cylindrical  Configuration. 


This  configuration  was  the  easiest  to  set  up  and  hence  was  done  first.  The  explosive 
chamber  and  explosive  within  it  were  provided  as  cylinders  and  therefore  could  be  left 
as  specified  above.  The  lengths  of  the  input  and  output  tunnels  were  retained  at  30  m, 
and  the  cross-sectional  areas  were  preserved  at  8  m  x  6  m  =  48  m2.  This  resulted  in  a  ra¬ 
dius  for  the  equivalent  cylindrical  input  and  output  tunnels  of  3.91  m.  The  expansion 
chamber  was  given  a  cylindrical  equivalent  in  the  same  manner,  retaining  the  length  of 
15  m  and  calculating  an  equivalent  radius  to  preserve  the  cross-sectional  area.  The  area 
was  70  m  x  8  m  =  560  m2,  thus  the  cylindrical  equivalent  radius  is  13.35  m.  The  final 
configuration  used  for  the  cylindrical  calculation  is  given  in  Figure  2. 


0.00 

I 


0.00 


80.20  [569  ZOOM)  95.20  (6/5  zones) 


2528  (180  zone*)  5020  (356  zones) 

J _ [  _  5.64  (40  zones) 

- -  1.41  (10  zones) 


1325  (95  zones) 

125.20  (888  zones) 

- -  —  3.91  (28  zones) 


Dimensions  In  meteis 


12  zones  beyond 
portal.  hcreaslng 
by  ratio  10% 


Configuration  for  CyBndricaJ  Test  Calculation 
Total  zones:  96x900  =>86.400 
Zone  size:  14.1  cm  except  in  region  beyond  portal 


Figure  2.  Configuration  and  zoning  used  for  cylindrical  calculation. 

1.2.2.  Two-Dimensional  Cartesian  Configuration. 

For  a  two-dimensional  Cartesian  calculation,  it  is  necessary  to  choose  a  "unit  height,"  or 
dimension  normal  to  the  two  dimensions  of  the  calculations!  mesh,  which  is  the  same 
for  every  part  of  the  calculation.  Because  the  input  and  output  tunnels  are  6  m  high,  we 
chose  6  m  as  the  unit  height.  The  dimensions  of  these  two  tunnels  thus  remained  file 
same,  and  those  of  the  detonation  and  expansion  chambers  were  adjusted  to  give  them 
equivalent  cross-sectional  area  using  the  6-m  unit  height.  The  lateral  dimension  of  the 
explosive  was  adjusted  in  the  same  way  to  retain  the  total  yield. 

The  cylindrical  explosive  chamber  has  a  radius  of  5.64  m,  making  its  cross-sectional  area 
99.93  m2.  Therefore,  for  a  6-m  unit  height  in  rectangular  coordinates,  its  width  must  be 
99.93  m2/6  m  =  16.66  m.  The  length  of  the  chamber  is  retained  at  50.20  m.  Similarly,  the 
expansion  chamber  has  a  cross-sectional  area  of  560  m2,  as  given  with  the  8-m  height. 
To  retain  the  length  and  volume  of  this  chamber  with  the  6-m  height,  we  have  to  in- 
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crease  its  width  from  70  m  to  93.33  m.  The  lateral  offset  between  the  input  and  output 
tunnels  was  retained  at  the  given  58  m,  so  the  additional  width  resulted  in  larger  al¬ 
coves  at  the  sides  of  the  expansion  chamber. 

The  dimensions  of  the  explosive  had  also  to  be  adjusted  so  that  the  proper  explosive 
yield  would  be  contained  in  the  Cartesian  problem.  The  volume  of  TNT,  as  calculated 
above,  is  158.52  m3.  Retaining  its  length  at  25.38  m  so  that  the  detonation  will  occur  at 
the  proper  rate,  and  using  the  6-m  unit  height,  the  width  of  the  explosive  slab  is  1.04  m. 
All  of  the  dimensions  used  in  this  configuration  are  shown  in  Figure  3.  Note  that  this  is 
a  plan  view,  so  that  discussions  of  width  or  lateral  dimension  are  shown  vertically  on 
the  paper,  while  length  is  shown  horizontally.  The  unit  height,  a  constant  6  m,  is  normal 


to  the  plane  of  the  paper. 


Figure  3.  Configuration  for  two-dimensional  Cartesian  problem  set-up. 

1.2.3.  Three-Dimensional  Configuration  (Cartesian). 

The  three-dimensional  problem  was  more  difficult  than  the  others  because  of  the  com¬ 
plexities  of  visualizing  and  displaying  the  three-dimensional  representation.  Because 
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everything  could  be  represented  exactly  as  it  was  provided,  this  phase  of  the  setup  was 
straightforward.  The  detonation  chamber,  input  tunnel,  and  explosive  were  all  cylindri¬ 
cal  and,  in  fact,  the  first  phase  of  the  problem  utilized  the  same  calculation  that  had 
been  generated  for  the  cylindrical  representation.  We  did  not  begin  the  three-dimen¬ 
sional  part  until  the  explosive  was  completely  detonated  and  the  shock  wave  had  pro¬ 
gressed  through  the  input  tube  and  was  beginning  to  move  into  the  expansion  chamber. 
At  this  point  the  three-dimensional  configuration  illustrated  in  Figure  4  was  set  up. 
Because  the  three-dimensional  code  allows  cylindrical  packages  in  the  Cartesian  mesh, 
representing  them  as  closely  as  possible  following  cell  boundaries,  no  modifications  to 
the  chamber  shapes  beyond  those  defined  by  the  "stair-step"  cell  boundaries  were  re¬ 
quired.  Figure  5  is  an  isometric  projection  of  the  information  provided  in  Figure  4,  and 
may  be  helpful  in  visualizing  the  configuration. 


Figure  4.  Configuration  drawing  for  three-dimensional  Cartesian  problem  set-up. 
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Figure  5.  Orthogonal  drawing  for  three-dimensional  configuration. 
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SECTION  2 

CALCULATION  SETUP 


2.1.  Cylindrical  Configuration. 


The  cylindrical  calculation  is  radially  symmetric  about  a  vertical*  axis  and  is  constructed 
uniformly  with  14.1  cm  zones  in  both  coordinate  directions.  For  the  charge  detonation, 
or  first  phase,  of  this  calculation,  the  fluid  region  extends  horizontally  to  a  radius  of  5.64 
m.  Two  additional  zones  of  perfectly  reflecting  solid  are  placed  to  the  right,  extending 
the  computational  domain  to  5.922  m.  Vertically,  the  domain  extends  to  50.196  m.  A 
cylindrical  region  with  a  radius  of  1.41  m  and  length  of  25.38  m  is  specified  as  TNT.  An 
explosive  mass  of  247  metric  tonnes  is  generated  which  contains  6.56x10^  ergs  of  am¬ 
bient  internal  energy.  The  detonation  of  this  cylinder  of  TNT  is  initiated  with  a  hemi¬ 
spherical  region  of  high-energy  gaseous  detonation  products.  This  region  is  19  cm  in 
radius  and  is  located  at  the  base  of  the  TNT  on  its  centerline.  The  initialization  hemi¬ 
sphere  radius  is  assigned  to  correspond  to  initialization  time  at  29  microseconds,  assum¬ 
ing  initiation  at  the  center  at  time  0  and  a  detonation  velocity  of  6.5  x  10^  cm/ s.  After  29 
microseconds,  detonation  proceeds  normally  through  the  TNT  according  to  well-estab¬ 
lished  bum  routines  in  the  hydrocode.  Figure  6  illustrates  the  configuration  of  the  initial 
cylindrical  detonation  calculation. 

The  data  representing  the  initial  conditions  as  described  was  processed  by  SHARC,  us¬ 
ing  second-order  alternating  direction  remapping,  a  multi-material  equation  of  state, 
and  a  high-explosive  bum  routine.  The  high-explosive  bum  was  complete  at  4  millisec¬ 
onds.  At  8  milliseconds,  the  shock  front  was  close  enough  to  the  top  boundary  of  the 
computational  domain  to  necessitate  remapping  the  calculation  to  a  larger  domain. 

The  larger  domain  consists  of  uniform  14.1  cm  zones  in  the  radial  direction  to  a  distance 
13.82  m.  In  the  vertical  direction,  the  zones  are  14.1  cm  uniformly  to  a  height  of  125.20 
m  and  expand  geometrically  beyond  that  point  by  10%  from  one  zone  to  the  next,  to  a 
final  height  of  128.20  m. 


*  The  "vertical"  and  "horizontal"  designations  here  and  in  Figures  6  and  7  are  reversed  from  those  shown 
in  Figures  1  and  2.  The  latter  were  drawn  to  correspond  to  physical  reality;  the  former  are  as  displayed  by 
the  hydrocode.  Because  there  is  no  treatment  of  gravity  in  these  calculations,  the  results  are  not  affected. 
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The  fluid  region  is  extended  into  three  additional  chambers.  The  first  extends  the  fluid 
region  30.03  m  to  a  height  of  80.22  m.  This  extension  is  constricted  to  a  width  of  3.948  m 
in  order  to  provide  a  chamber  of  comparable  cross-section  and  length  to  the  prescribed 
input  tube.  The  second  extends  the  fluid  region  14.94  m  to  a  height  of  95.17  m.  This 
chamber  expands  to  a  width  of  13.395  m  in  order  to  provide  the  appropriate  volume  for 
the  expansion  chamber. 

The  final  (output)  chamber  replicates  the  input  chamber  using  the  same  length  and 
width.  It  extends  the  fluid  region  to  a  height  of  125.208  m.  Beyond  this  position  are  ex¬ 
panded  zones  and  transmissive  boundaries  in  either  direction.  The  shock  may  expand 
without  restriction  in  this  region.  Figure  7  illustrates  the  configuration  for  the  remap  of 
the  cylindrical  calculation. 


Figure  6.  Configuration  for  initial  detonation  calculation  in  cylindrical  symmetry. 
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Figure  7.  Configuration  for  remap  of  cylindrical  calculation. 

2.2.  Two-Dimensional  Cartesian 

The  two-dimensional  Cartesian  calculation  consists  initially  of  a  domain  extending  hor¬ 
izontally  to  95.42  m  and  vertically*  to  35.36  m.  The  fluid  region  consists  of  two  horizon¬ 
tally-oriented  and  one  vertically-oriented  chambers.  Each  of  the  chambers  and  pas¬ 
sageways  is  assigned  a  unit  depth  of  6  m. 

The  chamber  in  which  the  detonation  is  initiated  is  16.64  m  wide  and  50.18  m  long.  It  is 
oriented  horizontally  with  its  centerline  vertically  offset  from  the  bottom  of  the  mesh  by 

*  Because  this  is  a  plan  view,  "vertical"  means  lateral  across  the  width  of  the  expansion  chamber,  just  as 
in  Figure  3. 
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17.68  m.  Specified  within  this  chamber,  along  its  centerline,  is  a  region  of  TNT  1.04  m 
wide  and  25.35  m  in  length.  This  rectangle  when  multiplied  by  the  6-m  unit  depth,  con¬ 
tains  a  volume  corresponding  to  247  metric  tonnes  of  TNT.  The  calculation  of  detona¬ 
tion  within  this  region  is  initiated  with  a  semi-circular  region  of  high  energy,  gaseous 
detonation  product  39  cm  in  diameter.  The  calculation  is  assigned  an  initial  time  of  59 
microseconds  to  correspond  to  this  diameter. 

Adjacent  to  the  detonation  chamber  is  a  constrictive  passageway  (the  input  chamber) 
30.0  m  long  and  8.06  m  in  width.  One  additional,  vertically-oriented,  rectangular 
chamber  is  located  to  the  right  of  the  passageway.  It  extends  the  fluid  region  14.95  m  to 
a  horizontal  distance  of  95.16  m  from  the  origin,  and  represents  part  of  the  expansion 
chamber.  Figure  8  illustrates  die  configuration  of  the  initial  Cartesian  detonation  calcu¬ 
lation.  The  zone  size  for  this  part  of  the  calculation  is  13  cm  in  both  coordinate  direc¬ 
tions. 

The  data  representing  the  initial  condition  as  described  was  processed  by  SHARC,  as  in 
the  previous  case.  The  high  explosive  bum  was  complete  at  approximately  4  millisec¬ 
onds.  At  18  milliseconds,  the  shock  front  was  close  enough  to  the  top  boundary  of  the 
computational  domain  to  necessitate  remapping  to  a  larger  domain. 

The  larger  domain  includes  an  extension  of  the  vertically-oriented  expansion  chamber 
to  93.34  m.  Additionally,  a  horizontally-oriented  output  chamber  30.0  m  by  8.0  m  is  lo¬ 
cated  adjacent  to  the  expansion  chamber.  Its  centerline  is  vertically  offset  75.66  m.  The 
zones  to  the  left  of  the  entrance  to  the  expansion  chamber  expand  geometrically  with  a 
10%  expansion  ratio  to  reduce  the  size  of  the  computational  files  and  processing  time. 
Figure  9  illustrates  the  configuration  of  the  final  remap  of  the  two-dimensional 
Cartesian  calculation. 
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Figure  8.  Configuration  for  initial  detonation  calculation  in  Cartesian  coordinates. 


Figure  9.  Configuration  for  remap  of  Cartesian  calculation. 


2.3.  Three-Dimensional  Cartesian  Calculation. 

The  three-dimensional  Cartesian  calculation  consists  initially  of  two  adjacent  half-cylin¬ 
ders  and  a  rectangular  prism  in  a  domain  extending  95.00  m  horizontally  and  35.36  m 
vertically.  The  cylinders  and  rectangular  prism  are  halved  in  order  to  exploit  symmetry 
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of  the  configuration  above  and  below  a  horizontal  mid-plane.  The  half-cylinders  were 
initialized  with  the  hydrodynamic  parameters  from  the  two-dimensional  cylindrical 
calculation  at  14  milliseconds.  The  two-dimensional  cylindrical  data  is  extrapolated  into 
the  three-dimensional  coordinate  system  by  using  a  rotation  transformation  algorithm. 
The  horizontal  (x-coordinate)  sizes  of  the  zones  are  geometrically  expanded  from  the 
entrance  of  the  rectangular  prism  leftward  to  the  origin  using  a  10%  expansion  ratio. 
Zones  within  the  rectangular  prism  are  15  cm  in  the  x-coordinate  dimension.  In  the  ver¬ 
tical  (z-coordinate)  and  page-normal  (y-coordinate)  the  zone  size  is  14.1  cm.  At  23  mil¬ 
liseconds,  the  shock  front  was  close  enough  to  the  top  boundary  of  the  computational 
domain  to  necessitate  remapping  the  calculation  to  a  larger  domain. 

The  next  larger  domain  was  extended  vertically  to  60.0  m.  The  zone  sizes  remained  un¬ 
changed.  At  34  milliseconds,  the  shock  front  position  necessitated  remap  to  a  still  larger 
domain.  The  next  remap  extended  the  domain  both  horizontally  and  vertically. 
Horizontal  (x-coordinate)  and  vertical  (y-coordinate)  zone  sizes  were  doubled,  to  30  cm. 
A  cylindrical  output  chamber  was  placed  adjacent  to  the  rectangular  prism  representing 
the  expansion  chamber  at  a  vertical  offset  of  64.0  m.  An  unrestricted  expansion  region 
was  included  beyond  the  output  chamber.  The  horizontal  zones  within  both  these  re¬ 
gions  expanded  geometrically  at  a  10%  ratio.  Figure  10  illustrates  the  final  configura¬ 
tion  of  the  three-dimensional  Cartesian  configuration. 
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SECTION  3 

CALCULATIONAL  RESULTS 

Figures  11  and  12  are  pressure  contour  plots  for  the  two-dimensional  cylindrical  and 
Cartesian  configurations  at  2  milliseconds.  A  separate  plot  for  the  three-dimensional 
configuration  has  not  been  included,  because  it  is  identical  to  the  cylindrical  configura¬ 
tion  at  this  early  time.  At  this  time,  the  shock  front  is  at  about  the  same  position  for  each 
case.  The  maximum  pressure  is  somewhat  lower  (about  20%)  in  the  Cartesian  chamber, 
probably  because  there  are  more  computational  zones  across  the  cylindrical  explosive, 
and  hence  the  pressure  is  not  relieved  as  quickly  as  in  the  Cartesian  case.  The  cylindrical 
result  is  probably  a  better  approximation  of  reality  during  the  explosive  detonation.  In 
both  cases,  reflected  shocks  are  proceeding  inward  from  the  chamber  walls. 

Figures  13  and  14  are  pressure  contour  plots  at  8  milliseconds.  These  plots  show  that  the 
2-D  Cartesian  shock  front  has  advanced  4  m  beyond  the  position  of  the  cylindrical 
shock.  The  Cartesian  shock  front  is  already  reflecting  off  the  end  wall  of  the  detonation 
chamber.  The  shock  structures  differ  significantly  because  of  the  geometry.  In  the  cylin¬ 
drical  configuration,  a  Mach  stem  is  evident  on  the  axis,  whereas  no  such  structure 
shows  in  the  Cartesian  contour  plot.  The  maximum  pressure  is  about  a  factor  of  2.5 
stronger  in  the  cylindrical  chamber.  These  effects  are  almost  surely  a  result  of  cylindri¬ 
cal  convergence  on  the  axis  of  waves  reflected  at  the  chamber  walls.  The  Cartesian  con¬ 
figuration  does  not  generate  geometric  convergence. 

Figures  15, 16,  and  17  are  pressure  contour  plots  at  16  milliseconds.  Because  the  2-D 
cylindrical  configuration  was  mapped  into  the  3-D  mesh  at  14  milliseconds,  the  3-D 
Cartesian  is  a  distinct  calculation  at  16  milliseconds  and  is  included  for  comparison.  The 
contours  shown  are  on  the  horizontal  reflecting  plane  at  the  center  of  the  configuration. 
At  16  milliseconds,  the  shock  wave  is  entering  the  expansion  chamber  in  each  case.  The 
2-D  Cartesian  case  has  an  unrealistically  long  alcove  at  the  end  of  the  expansion  cham¬ 
ber  (bottom  right  of  the  plot),  made  necessary  by  the  scheme  used  to  maintain  the  same 
expansion  chamber  volume  in  all  three  cases.  This  means  that  the  end  wall  reflection  oc¬ 
curs  much  later  in  this  case  than  in  the  3-D  case.  Note  that,  at  this  time,  the  maximum 
pressure  occurs  in  the  2-D  Cartesian  case. 

Figures  18, 19,  and  20  are  pressure  contour  plots  at  similar,  but  not  exactly  the  same, 
times.  The  Cartesian  plots  are  taken  at  26  ms  for  the  3-D  and  at  25  milliseconds  for  the  2- 
D.  The  shock  structures  are  similar  at  these  times.  For  the  Cartesian  cases,  the  shock 
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front  is  a  diagonal  across  the  expansion  chamber,  arriving  earlier  at  a  position  close  to 
the  far  wall  than  at  a  corresponding  position  on  the  near  wall.  The  cylindrical  case  can¬ 
not  be  easily  compared  to  the  others  at  these  times,  because  the  shock,  having  a  much 
shorter  distance  to  travel  through  the  expansion  chamber,  has  gone  on  into  the  output 
tube  and  is  near  the  tunnel  portal. 

Figures  21  and  22  are  pressure  contour  plots  showing  shock  front  development  at  36 
milliseconds.  The  2-D  Cartesian  case  has  a  shock  front  a  bit  further  into  the  output 
chamber  than  does  the  3-D  case.  This  is  because  the  3-D  expansion  chamber  is  taller  (8 
m  as  opposed  to  the  6-m  unit  height  of  the  2-D  configuration),  and  hence  the  shock  front 
has  more  room  in  which  to  expand.  This  will  allow  the  shock-front  pressure  to  decay 
more  rapidly  and  hence  the  front  will  not  travel  as  fast.  This  volume  effect  is  eventually 
equalized,  but  not  until  the  wave  has  completely  filled  the  alcove  at  the  end  of  the 
chamber.  We  do  not  show  a  cylindrical  contour  at  this  time  because  for  this  case,  the 
shock  has  exited  the  output  tunnel. 

The  next  group  of  figures  consists  of  comparisons  of  overpressure  records  obtained 
from  the  three  calculations  at  corresponding  points  in  the  interior  and  at  the  portal. 
Figure  23  provides  a  comparison  of  the  overpressure  signal  at  the  detonation  chamber 
exit.  Only  the  2-D  cylindrical  and  the  2-D  Cartesian  are  compared  because  the  3-D 
Cartesian  is  identical  to  the  2-D  cylindrical  at  this  location.  The  signal  beyond  the  initial 
shock  arrival  peaks  is  roughly  comparable  with  respect  to  amplitude  and  general  wave¬ 
form.  Differences  in  timing  and  amplitude  of  the  various  peaks  occur  because  of  differ¬ 
ences  in  the  explosive  configuration  and  in  the  distances  and  orientations  of  various  re¬ 
flecting  surfaces  relative  to  the  measurement  location. 

Figure  24  compares  overpressure  waveforms  at  locations  in  the  center  of  the  input  tube. 
As  expected,  the  arrival  times  are  nearly  identical.  The  peak  overpressure  is  somewhat 
lower  in  the  Cartesian  case,  but  the  generation  of  impulse,  as  shown  in  Figure  25,  is  ap¬ 
proximately  equivalent.  Similar  results  at  the  input  tube/ expansion  chamber  junction 
are  shown  in  Figures  26  and  27. 

Figure  28  compares  the  signals  obtained  near  the  center  of  the  expansion  chamber.  The 
three-dimensional  result  is  included  for  this  and  the  following  comparisons  because  at 
this  location,  the  record  is  different  from  that  of  the  cylindrical  calculation.  The  cylindri¬ 
cal  record  exhibits  an  earlier  arrival  time  because  the  path  over  which  the  shock  wave 
travels  is  shorter.  It  is  also  greater  in  amplitude.  The  Cartesian  cases  appear  to  have 
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similar  magnitudes  for  the  initial  5  milliseconds  after  shock  arrival.  The  late-time  signal 
is  stronger  in  the  3-D  case,  possibly  because  of  the  larger  alcove  volume  of  the  expan¬ 
sion  chamber  in  the  2-D  configuration. 

Figure  29  compares  signals  from  the  center  of  the  output  tube.  Again,  the  cylindrical 
signal  arrives  much  earlier  and  has  higher  peaks.  The  Cartesian  cases  are  in  closer 
agreement,  although  there  are  noticeable  differences  in  waveform  details.  The  same 
may  be  observed  at  the  exterior  opening  of  the  output  tube,  as  shown  in  Figure  30. 
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Figure  11.  Pressure  contour  plot  for  two-dimensional  cylindrical  detonation  calculation 
at  2  ms. 
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Figure  12.  Pressure  contour  plot  for  two-dimensional  Cartesian  detonation  calculation 
at  2  ms. 
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Figure  13.  Pressure  contour  plot  for  two-dimensional  cylindrical  calculation  at  8  ms. 
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Figure  14.  Pressure  contour  plot  for  two-dimensional  Cartesian  calculation  at  8  ms. 
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Figure  15.  Pressure  contour  plot  for  three-dimensional  Cartesian  calculation  at  16  ms. 
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Figure  16.  Pressure  contour  plot  for  two-dimensional  Cartesian  calculation  at  15  ms. 
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Figure  17.  Pressure  contour  plot  for  two-dimensional  cylindrical  calculation  at  16  ms. 
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Figure  18.  Pressure  contour  plot  for  three-dimensional  Cartesian  calculation  at  26  ms. 
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Figure  19.  Pressure  contour  plot  for  two-dimensional  Cartesian  calculation  at  25  ms. 
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Figure  21.  Pressure  contour  plot  for  three-dimensional  Cartesian  calculation  at  36  ms 


Z 


ALTITUDE  U 


PRESSURE 


RANCE  U 

CHARGE  DETONATION  VENTED  CHAMBER  (CARTESIAN)  RGE  19  OCT  1.15 
TIME  36.000  MSEC  CYCLE  1763.  PROBLEM  281.0011 

Figure  22.  Pressure  contour  plot  for  two-dimensional  Cartesian  calculation  at  36  ms. 
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Figure  23.  Comparison  of  overpressure  waveforms  at  entrance  to  input  tube  from  two- 
dimensional  (a)  cylindrical  and  (b)  Cartesian  calculations. 
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Figure  24.  Comparison  of  overpressure  waveforms  at  center  of  input  tube  from  two-di¬ 
mensional  (a)  cylindrical  and  (b)  Cartesian  calculations. 
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Figure  25.  Comparison  of  overpressure  impulse  records  at  center  of  input  tube  from 
two-dimensional  (a)  cylindrical  and  (b)  Cartesian  calculations. 
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Figure  26.  Comparison  of  overpressure  waveforms  at  entrance  to  expansion  chamber 
from  two-dimensional  (a)  cylindrical  and  (b)  Cartesian  calculations. 
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igure  27.  Comparison  of  overpressure  impulse  records  at  entrance  to  expansion  cham¬ 
ber  from  two-dimensional  (a)  cylindrical  and  (b)  Cartesian  calculations. 
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Figure  28.  Comparison  of  overpressure  waveforms  at  center  of  expansion  chamber 
from  two-dimensional  (a)  cylindrical  and  (b)  Cartesian  calculations,  and 
from  (c)  three-dimensional  Cartesian  calculation. 


OVER  PRESSURE 


rwAPOr  DETONATION  VENTED  CHAMBER  (CYLINDRICAL)  RGE 
282J300I-STA-52 


OVERPRESSURE 


CHARGE  DETONATION  VENTED  CHAMBER  (CARTESIAN)  RGE 
281X010.STAJ8 


O  JOOjO 


OVERPRESSURE 


S 

a 

i 

or 

\ 

H 

Z 

a 

w 

« 

CO 

w 

w 

c: 

cu 


400.0  - 

3000- 

2000- 

1000- 

oo- 


.  I  1,1 

Si)  loo 


i 

15i)  200  2SJ3  300 


TIMS  (SECONDS  X  10" 


3) 


400 


- - r 

45  J)  500 


CHARGE  DETONATION  VENT  sD  CHAMBER  (WES-CRAY-3D)  RGE 
283.0003.STA.73 


Figure  29.  Comparison  of  overpressure  waveforms  at  center  of  output  tube  from  two 

dimensional  (a)  cylindrical  and  (b)  Cartesian  calculations,  and  from  (c)  three- 
dimensional  Cartesian  calculation. 
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Figure  30.  Comparison  of  overpressure  waveforms  at  exterior  end  of  output  tube  from 
two-dimensional  (a)  cylindrical  and  (b)  Cartesian  calculations,  and  from  (c) 
three-dimensional  Cartesian  calculation. 
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SECTION  4 


COMPARISON  OF  APPROACHES 

4.1.  Arrival  Time,  Peak  Overpressure,  and  Peak  Dynamic  Pressure  Comparisons 

Figures  31, 32,  and  33  are  comparisons  of  arrival  times,  peak  overpressures  and  peak 
dynamic  pressures  at  selected  points  within  the  tunnel  complex.  For  the  independent 
variable,  we  use  "cumulative  distance,"  which  is  defined  as  the  distance  from  the  center 
back  of  the  detonation  chamber  (the  detonation  point)  along  an  orthogonal  path  to  the 
point  in  question  in  the  Cartesian  cases.  The  path  follows  the  centerline  of  the  detona¬ 
tion  chamber  and  input  tube  to  the  center  of  the  expansion  chamber  at  the  input  end.  It 
then  makes  a  right-angle  turn  and  continues  through  the  expansion  chamber  for  the 
length  of  the  input  tube/output  tube  offset  (58  m).  A  second  right  angle  turn  directs  the 
path  through  the  output  tube  to  the  portal.  Key  points  on  the  cumulative  distance  path 
are  given  below; 


Key  Point  on  Cumulative 

Cumulative  Distance  to  This 

Distance  Path 

Point 

Detonation  Chamber /Input 

50.2  m 

Tube  Interface 

Input  Tube/Expansion 

80.2  m 

Chamber  Interface 

Center  of  Expansion  Chamber 

87.7  m 

at  Input  Tube  End 

Center  of  Expansion  Chamber 

116.7  m 

Center  of  Expansion  Chamber 

145.7 

at  Output  Tube  End 

Expansion  Chamber/Output 

153.2  m 

Tube  Interface 

Portal 

183.2  m 
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For  the  cylindrical  results  given  in  the  plots,  the  independent  variable  has  been  adjusted 
so  that  key  points  in  the  cylindrical  configuration  are  plotted  at  the  corresponding 
points  listed  above.  The  cylindrical  results  are  identical  to  those  for  the  3-D  Cartesian 
results  in  the  detonation  chamber  and  input  tube;  thus  only  the  2-D  and  3-D  Cartesian 
results  are  plotted  in  this  region. 

The  Cartesian  arrival  times.  Figure  31,  are  in  good  agreement  throughout  the  calcula¬ 
tions.  The  3-D  shock  is  slightly  later  in  the  expansion  chamber  and  beyond,  probably 
because  the  higher  overhead  (8  m  compared  to  the  6-m  unit  height  in  the  2-D  case)  al¬ 
lows  more  local  expansion  volume  as  the  shock  travels  through  the  expansion  chamber. 
The  cylindrical  case  arrival  time  is  much  earlier,  as  expected,  because  of  the  shorter  dis¬ 
tance  of  travel  to  corresponding  points  in  the  complex. 

The  peak  overpressure  plot.  Figures  32,  shows  similar  values  for  the  two  Cartesian 
cases;  however  the  3-D  values  are  slightly  higher  in  all  cases.  This  might  be  because  the 
explosive  detonation  was  better  characterized  by  the  axisymmetric  representation.  We 
expected  the  2-D  result  to  be  almost  the  same  as  the  3-D  during  transition  of  the  expan¬ 
sion  chamber ,  but  it  never  quite  achieves  the  same  amplitude.  The  cylindrical  result  is 
higher  by  nearly  an  order  of  magnitude  because  there  are  no  right-angle  turns  to  im¬ 
pede  the  flow. 

This  problem  caused  by  unimpeded  axial  flow  in  the  cylindrical  case  is  particularly  no¬ 
ticeable  in  Figure  33,  in  which  peak  dynamic  pressures  are  plotted.  For  this  plot,  the  dy¬ 
namic  pressure  component  shown  for  comparison  was  that  in  the  primary  direction  of 
flow.  For  the  cylindrical  case,  this  was  along  the  axial  direction.  For  the  Cartesian  con¬ 
figurations,  we  used  the  horizontal  (in  the  previously  described  calculational  represen¬ 
tation)  component  within  the  detonation  chamber  and  input  and  output  tubes,  the  ver¬ 
tical  component  in  the  expansion  chamber.  Because  there  is  crossflow  in  all  the  cham¬ 
bers,  this  may  not  have  been  the  best  choice,  and  may  be  responsible  for  the  erratic  be¬ 
havior  of  the  curves. 
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Cumulative  Distance  (m) 


Figure  31.  Arrival  time  comparison  for  the  three  configurations  along  an  orthogonal 
central  path  from  detonation  point  to  portal. 
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Figure  32.  Peak  overpressure  comparison  for  the  three  configurations  along  an  or¬ 
thogonal  central  path  from  detonation  point  to  portal. 
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Figure  33.  Peak  dynamic  pressure  (in  primary  flow  direction)  for  tire  three  configura¬ 
tions  along  an  orthogonal  central  path  from  detonation  point  to  portal. 


4.2.  Discussion 

The  three  geometric  representations  were  shown  approximately  in  Figure  1,  and  in 
more  detail  in  the  subsequent  figures.  The  three-dimensional  case  has  a  cross-section 
reduction  between  the  detonation  chamber  and  connecting  tunnel  (or  input  tube),  a  sig¬ 
nificant  expansion  in  cross  section  and  a  right-angle  turn  upon  transition  to  the  expan¬ 
sion  chamber,  a  travel  distance  of  58  m  before  another  right-angle  turn  and  a  cross-sec¬ 
tion  reduction  to  the  output  tube.  The  two-dimensional  cylindrical  approximation  in¬ 
cludes  the  appropriate  cross-sectional  area  changes  but  has  no  changes  of  direction.  The 
total  length  of  the  expansion  chamber,  as  traversed  by  the  shock  wave,  is  considerably 
shorter  in  the  cylindrical  than  in  the  other  cases.  Therefore  we  expected,  and  observed, 
reduced  arrival  times  in  the  cylindrical  case.  Because  there  were  no  changes  in  direction, 
the  dynamic  pressures  were  significantly  higher  than  in  either  of  the  Cartesian  cases. 
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This  lack  of  direction  changes  was  also  a  factor  contributing  to  the  reduced  arrival 
times. 

In  the  2-D  Cartesian  case,  the  right-angle  turns  and  the  distance  between  turns  was  ac¬ 
curately  represented.  The  detonation  chamber,  originally  cylindrical,  was  modified  to 
equalize  its  volume.  Also,  the  expansion  chamber  shape  was  modified.  The  cross-sec¬ 
tional  area  of  the  expansion  chamber,  for  a  plane  normal  to  the  major  flow  direction, 
was  somewhat  less  than  in  the  true,  3-D  geometry.  The  resulting  difference  in  volume 
was  accounted  for  by  extending  the  expansion  chamber  equally  on  each  side  beyond  the 
input  and  output  tube  connections.  These  changes  account  for  the  differences  noted  in 
these  regions.  Because  the  volumes,  distances  and  angles  in  the  2-D  and  3-D  Cartesian 
cases  are  nearly  equal,  the  results  from  the  two  calculations  at  the  portal  are  reasonably 
close. 


42 


SECTION  5 

CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  study,  we  set  up  the  same  basic  situation  in  three  different  calculational  configu¬ 
rations,  performed  the  calculations,  and  then  compared  results  to  determine  whether  a 
two-dimensional  approximation  is  valid  to  represent  a  three-dimensional  problem,  and 
if  so,  whether  a  cylindrical  or  Cartesian  representation  most  nearly  reproduces  the  true 
three-dimensional  result.  The  answer  is,  and  this  might  have  been  expected,  that  any 
time  a  change  is  made  to  the  configuration,  this  change  will  be  reflected  to  some  degree 
in  the  emerging  waveform.  Ideal  surfaces  reflect  shocks,  and  if  these  surfaces  are  at  dif¬ 
ferent  distances  from  the  source,  or  are  oriented  differently,  then  the  reflected  waves,  as 
they  contribute  to  the  form  of  an  out-flowing  wave,  will  arrive  at  different  times  and 
have  different  characteristics.  This  cannot  be  avoided. 

In  many  cases,  however,  the  exact  waveform  to  be  observed  at  some  point  within  or  just 
outside  of  a  tunnel  complex  is  not  as  important  as  the  peak  overpressure,  peak  dynamic 
pressure,  or  their  respective  impulse  values.  In  these  cases,  two-dimensional  represen¬ 
tations,  if  they  are  carefully  done,  can  provide  adequate  predictions  of  test  cases  or  ade¬ 
quate  simulations  for  detonation  scenarios  of  interest. 

A  three-dimensional  simulation  gives,  of  course,  the  closest  correspondence  to  the  de¬ 
sired  situation.  This  must  be  balanced  by  the  fact  that  the  third  dimension  may  require 
one  to  two  orders  of  magnitude  more  computing  time,  because  the  computer  is  doing 
the  equivalent  of  a  number  of  two-dimensional  calculations  equal  to  the  number  of  cells 
in  the  third  direction.  A  large,  fast,  mainframe  computer  is  usually  required  to  complete 
a  three-dimensional  calculation,  and  these  are  not  always  readily  available.  A  compro¬ 
mise  is  often  reached  by  using  a  coarser  mesh  size,  which  may  or  may  not  be  appropri¬ 
ate,  depending  on  requirements  of  the  calculational  results. 

Based  on  the  results  of  this  study,  it  appears  that  a  two-dimensional  Cartesian  represen¬ 
tation  is  a  reasonable  choice  for  most  tunnel  configurations  which  include  comers 
around  which  the  shock  wave  must  flow.  Such  a  configuration  provides  a  reasonably 
accurate  simulation  of  the  actual  interior,  although  problems  may  arise,  for  example,  if 
the  explosive  charge  is  fairly  small  so  that  stretching  it  over  the  whole  unit  height  makes 
it  ridiculously  narrow.  It  has  been  our  experience  that  a  charge  which  is  less  than  about 
five  computational  cells  wide  will  not  bum  properly,  because  of  the  pressure  release  at 
the  edges. 
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The  cylindrical  configuration  is  not  really  appropriate  for  configurations  with  comers 
like  the  one  of  interest  here.  The  path  that  the  shock  wave  follows  in  the  cylindrical  case 
is  a  straight  shot  along  the  centerline;  it  is  not  degraded  by  bouncing  off  of  walls  as  it 
goes  around  comers.  The  path  lengths  are  not  correct,  so  the  times-of-arrival  are  not 
correct.  This  configuration  should  provide  a  "worst-case"  answer,  but  it  may  be  too 
conservative  in  many  instances. 

While  analyzing  the  results  of  the  cylindrical  calculation,  we  thought  of  a  modification 
that  might  be  incorporated  into  the  configuration  to  provide  a  better  simulation.  This 
consists  of  placing  a  "blocker"  in  the  expansion  chamber  on  the  centerline  to  break  up 
the  flow,  as  shown  in  Figure  34.  The  outer  diameter  of  the  expansion  chamber  would 
have  to  be  modified  slightly  to  maintain  the  appropriate  expansion  chamber  volume, 
and  the  size  of  the  blocker  could  be  modified  to  control  total  path  length,  and  hence  ar¬ 
rival  time  of  the  shock  at  the  exit  tunnel.  Care  would  need  to  be  taken  not  to  choke  the 
flow  by  making  the  entrance  and  exit  to  the  expansion  chamber  smaller  in  the  cross  sec¬ 
tion  than  the  entrance  and  exit  tubes  themselves. 


1 

i 

i  blocker 
i 

Figure  34.  Cylindrical  expansion  chamber  configuration  with  "blocker". 

The  Cartesian  two-dimensional  calculation  could  also  have  been  set  up  in  a  slightly  dif¬ 
ferent  manner  which,  in  retrospect,  might  have  been  more  realistic.  Referring  to  Figure 
3,  recall  that  the  "length"  of  the  expansion  chamber,  shown  as  the  horizontal  dimension 
in  the  figure,  was  actually  the  shortest  dimension  of  the  chamber.  This  was  held  con¬ 
stant,  while  the  "width"  (shown  vertically)  was  adjusted  to  maintain  the  appropriate 
volume.  This  approach  required  the  addition  of  alcoves  at  the  sides  of  the  expansion 
chamber  whose  ends  were  farther  from  the  entrance  and  exit  tunnels  than  those  in  the 
three-dimensional  configuration,  and  hence  modified  the  timing  of  waves  reflected 
from  the  side  walls.  Another  approach  would  have  been  to  maintain  the  width  constant, 
and  to  make  an  adjustment  in  the  length  of  the  expansion  chamber.  This  modification  is 
illustrated  in  Figure  35,  which  shows  the  changes  to  the  expansion  chamber  shape  in 


44 


heavy  dashed  lines.  The  "width"  (vertical)  of  the  expansion  chamber  is  maintained  at  70 
m,  while  the  'length"  (horizontal)  is  increased  to  20  m  to  accommodate  the  6-m  unit 
height.  The  volume  is  the  same  in  both  cases. 


Figure  35.  Expansion  chamber  in  two-dimensional  Cartesian  system  illustrating  a 
possible  modification  to  the  chamber  shape. 
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